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Abstract 

We present a two-parameter family of recursion formulas for scalar field the- 
ory. The first parameter is the dimension (D). The second parameter (£) 
allows one to continuously extrapolate between Wilson's approximate recur- 
sion formula and the recursion formula of Dyson's hierarchical model. We 
show numerically that at fixed D, the critical exponent 7 depends continu- 
ously on (. We suggest the use of the (— independence as a guide to construct 
improved recursion formulas. 
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The renormalization group |],|2| method has shed light on the related questions of second 
order phase transitions and the infinite cut-off limit of field theories. The practical appli- 
cation of the method usually requires approximations. In the case of scalar field theory, a 
particularly simple approximate renormalization group transformation [Jj is Wilson's "ap- 
proximate recursion formula" (WARF in the following). The WARF is a simple integral 
equation with only one variable and has a free parameter which can be adjusted to ap- 
proximate a D-dimensional theory. The WARF can be handled very easily using numerical 
methods or various perturbative expansions. The numerical value of the critical exponent 
7 obtained with the the WARF is 1.218 for D = 3. The WARF can also be used to 
study non-perturbatively the cut-off dependence of the renormalized quantities. In the case 
D = 4, this can be used to set triviality bounds on the mass of a scalar particle. 

The derivation of the WARF is a masterpiece of quantum calculation. Unfortunately, 
it is not based on an expansion in a small parameter and there no obvious way to restore 
order by order the details erased by the approximation. Clearly, one needs an organizing 
principle to improve the WARF. A group theoretical approach |E| was proposed in the case 
of theories with quadratic interactions but this approach fails to control the proliferation of 
non-local terms in the case of quartic interactions. At the end of this letter, a new method 
of improvement will be suggested. 

A renormalization group transformation "very closely connected" || to the WARF holds 
exactly for Dyson's hierarchical model In the following, we use DHMRF as short for 

"Dyson's hierarchical model recursion formula" . The DHMRF is also an integral equation 
with one variable, and it also has a free parameter expressible in terms of D. For definiteness, 
the WARF and DHMRF are given below by Eqs. (7-8) with ( = 1 and 1/D respectively. 
Roughly speaking [7[], the WARF does in one step what the DHMRF does in D steps. 
However, this is not an exact statement because |9| the value of 7 in D = 3 is 1.300. 
Consequently, the two models have different physical properties. 

Before going further, we would like to emphasize that the difference between the two 
values of 7 quoted above is significantly larger than the errors involved in the numerical 
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calculations. We have repeated Wilson's calculation with smaller integration steps, cutting 
the integral at larger values and using different criteria to determine criticality. We found 
that these changes affected 7 by less than 0.002. We also used Wilson's numerical integration 



method for the DHRMF and found a value of 1.301 for 7. We confirmed [JTUJ] this result with 
errors smaller than 0.003 using the first 800 coefficients of the high-temperature expansion. 
Clearly, the difference between 1.22 and 1.30 is more than 25 times larger than the numerical 
errors involved in each calculation. 

The WARF and the DHMRF can be seen as two approximate versions of the (much more 
complicated) renormalization group transformation for a scalar lattice model with nearest 
neighbor iterations. The WARF integrates 2 D field variables at a time keeping their sum 
constant while the DHMRF does the same thing but with only two field variables. The fact 
that different values of 7 are obtained with the two approximate methods sets a limit on 
the accuracy of the approximation. For comparison, the universal value of 7 in D = 3 for 
nearest neighbor models is approximately 1.25. 

Ideally, one would like to construct a renormalization procedure for the nearest neighbor 
models corresponding to the integration of 

a = 2< D (1) 

field variables at a time and have all the physical quantities independent of £■ The quantity 
2^ is a scale factor which would play a role similar to the arbitrary scale parameter /1 used 



m some version 



|TT|j of the Callan-Symanzik equations and where one obtains homogeneous 
differential equations of the form ^ (physical quantity) =0. In this letter, we make a first 
step in a similar direction, by introducing a recursion formula where ( is arbitrary and which 
interpolates continuously between the WARF (( = 1) and the DHMRF (( = l/D). 

This recursion formula can be constructed using Dyson's hierarchical model (DHM) as 
an explicit realization and then extending the results for arbitrary In other words, this 
paragraph should first be read with ( taking the fixed value l/D and considered as well- 
known results written in a slightly unusual way. We call the sum of the fields in a cube 



containing a 1 sites fa, with I = 1,2.. corresponding to successive renormalization group 
transformations. We recall that a is defined in Eq.(l) and takes the value 2 for DHM. At 



criticality, one has the scaling law [12 



< (fa) 2 > oc (ab) 21 (2) 

where 

6 = 2"^. (3) 

Calling Pi(x)dx the probability for (fa/(ab) 1 ) to take a value between x and x + dx, we see 
that at criticality, 

pi = J dxPi(x)x 2 (4) 

tends to a constant for large /. On the other hand, in the high-temperature phase, 

< (fa) 2 > oc (a) 1 (5) 

and pi oc (ab 2 )~ l . Pi(x) is the main object studied with the renormalization group method. 
Using the parametrization, 

Pi{x) = Ke- (1/a)Q,(bx \ (6) 
one can check that the recursion formula 

Q l+ i(x) = -a H 1 -^} (7) 

with 

I t ( x ) = f°° dye- y2 -^ x+y) -^ Ql ^ y) (8) 

J — oo 

is equivalent to the DHMRF. The equivalence with the formulation usually found in the 
literature || is made clear by writing the recursion formula for the quantity ipi( x ) — h-i(bx). 
It is also clear that Eqs.(7-8) can be used for arbitrary value of (. In the case, ( = 1 



one recovers immediately the WARF. In conclusion, Eqs. (7-8) can be used to interpolate 
continuously between the WARF and the DHMRF. 

Note that b in Eq.(3) can be seen as the scaling factor of a massless Gaussian field under a 
change of scale 2 < \ This remark can be understood better by noting that the continuous (i.e. 
unregularized) version [13] of the massless Gaussian DHM is invariant under certain scale 



transformations. However, it is a discrete scale invariance. It seems plausible that discrete 
scale invariance allows log-periodic corrections 0] to the scaling laws. Such a corrections 



were indeed observed [|n| very clearly in the case of DHM. Or goal is to find a formulation 
where the discrete scale invariance (and the unphysical features associated with it) would 
be replaced by a continuous scale invariance. 

The critical behavior associated with the general recursion formula (7-8) can be studied 
with the usual methods. One starts with an initial function Qo(x) = rx 2 + gx 4 and for fixed 
g determines the critical value r c by observing the transition in the behavior of pi defined 
above. Near this critical value, one obtains the linearized expression 

Qi(x) = Q c {x) + {r-r c )\ l R c {x) (9) 

where A is the largest eigenvalue of the renormalization group transformation. From this, 
one finds the critical exponent 

The numerical values for various values of £ are displayed in Fig. 1. The results indicate 
that 7 is a continuous function of (. As ( becomes smaller, a larger number of iterations is 
necessary in order to obtain the critical value of r with an acceptable accuracy. For practical 
purposes, the number of iterations is of the order of 15/ The limit — is thus difficult to 
reach computationally, however we have found no indication of a drastic change of behavior 
(e.g., a sudden drop to 1) in this limit. 

We want to modify Eqs. (7-8) in order to get (- independent physical quantities. In the 
following, we focus the discussion on the linear behavior given by Eq. (9). First, the 



(■-independence of 7 requires that A(() = (A(l))^. In addition, the renormalized coupling 
constants should also be (-independent. Defining them with the procedure of Ref. [0 (which 
can be extended straightforwardly for arbitrary a and b), we see that the functions Q c and 
R c should also be (-independent. This requirement gives useful information concerning the 
corrections that need to be made to Eqs. (7-8). For instance if Q c is a (-independent fixed 
point of such a corrected formula, one can set ( = (0 + 5 in the known part of the recursion 
formula. Expanding in 5, one can construct, order by order in 5, "counterterms" that cancel 
the (-dependence. For instance at first order in S we need to add corrections which have a 
form similar to the r.h.s of (8) but with insertions of Qi and its first derivative with easily 
calculable coefficients. The effects of these corrections will be investigated using numerical 
methods. We expect that this procedure can be used to systematically improve the WARF 
and applied to realistic calculations of the critical exponents and the triviality bounds. 
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Figure Captions 

Fig. 1: The critical exponent 7 as a function of £ for D = 3. 



8 



1 — 1 — ' — 1 — I — 1 — ' — ' — I — ' — 1 — 1 — I — r 



— • 



CO 



_L 



00 



to 



d 



00 

d 

d 

to 
d 

d 

d 

CO 

d 

c\? 
d 



